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Carrier dynamics at energy level anticrossings in biased semiconductor superlattices, was studied 
in the time domain by solving the time-dependent Schrodinger equation. The resonant nature 
of interminiband Rabi oscillations has been explicitly demonstrated to arise from interference of 
intrawell and Bloch oscillations. We also report a simulation of direct Rabi oscillations across three 
minibands, in the high field regime, due to interaction between three strongly coupled minibands. 

PACS numbers: 73.23.-b, 73.21.Cd, 78.20.Bh, 78.30. Fs 



I. INTRODUCTION 

The rapid progress of solid-state electronics in recent 
years has drawn attention to semiconductor (SC) super- 
lattices (SL) as a useful source of coherent electrons. The 
proposed applications include microwave radiation (the 
so-called Bloch oscillator systems [H, 0, matter- wave 
interferometry Q and operational qubits for quantum 
computing [5|. A better understanding of fundamental 
types of carrier behaviour in biased superlattices is es- 
sential for making further progress in this field. 

Under bias, semiconductor superlattices demonstrate 
some remarkable quantum transport effects, such as res- 
onant Zener tunneling (RZT) and interminiband Rabi os- 
cillations [1, 0i Hj 0] • The phenomenon of interminiband 
Rabi oscillations in SC SL is also known in the literature 
as excitonic Rabi oscillations, Rabi flopping, periodic 
population swapping, field-induced delocalization, oscil- 
latory dipole, interwell oscillations and Bloch-Zener oscil- 
lations. Previous investigations [10,, [ll|, [l^] have covered 
many aspects of Rabi oscillations in SC systems since 
their first experimental observation in the early 1990's, 
by the laser pump-and-probe technique [l^. In experi- 
ment, Rabi oscillations occur as oscillating charge den- 
sity dipoles that have been observed both at low [ij] and 
room [l3| temperatures. A typical system demonstrat- 
ing Rabi oscillations are SC quantum dots under pulsed 
resonant excitation. 

While a two-miniband model works well for shallow 
superlattices (e.g. optical potentials), stronger potentials 
generally require a more elaborate approach. In the past 
few years some authors have employed more powerful 
calculational techniques, usually in the context of driven 
vertical transport, basing their work on resonant states 
and resonant Wannier-Stark functions of a system [l^, 
[T6j : however, the physical mechanism underlying Rabi 
oscillations has not been sufficiently elucidated. 

We avoid many simplifying yet restrictive model condi- 
tions, by directly solving the time-dependent Schrodinger 
equation along with transparent boundary conditions 
(TrBC). This enables us to go beyond the common two- 
band approximation in the high-field regime, thereby ob- 
taining a more reliable description of carrier dynamics, 
in particular Rabi oscillations, and allowing us to model 



new phenomena in quantum transport. This paper is 
organized as follows. Section [ll] provides details of the 
physical model used and its numerical implementation; 
the simulation results are presented in sections IIIIl IIVI 
and |Vl Section IIIIl deals with the occurrence and struc- 
ture of resonances; section ITVl discusses their nature, and 
self- interference of a wavepacket. Finally section |V] de- 
scribes the carrier dynamics at a resonance across three 
minibands as revealed in our simulations. 



II. PHYSICAL MODEL 

We consider a layered GaAs/Gai_j;Al3;As heterostruc- 
ture. The longitudinal motion of a wavepacket 5'(x,t) 
representing a single electron in a zero-temperature bi- 
ased superlattice with potential Vsl{x) under constant 
uniform bias F = —eE is described by a solution of the 
time-dependent Schrodinger equation 

- *(^' ^) + i^SLix) + xf) ^{x, t) 

= E^ix,t) (1) 

which we solved in the time domain. 

Experimentally the excitonic carrier population in 
the conduction miniband is created by ultrashort laser 
pulses [17]. This work deals only with conduction mini- 
band electrons, due to the fact that holes with their large 
effective mass are well-localized and do not demonstrate 
field-dependent absorption spectra [isj . 

A. Transparent boundary conditions 

Unlike previous studies based on a similar ap- 
proach [11,113,1211, we used TrBC [H,!!!,!!! that were 
recently derived for the Schrodinger equation in ID. We 
followed Moyer [25j] who used Crank-Nicholson to ad- 
vance the time, and the Numerov method for the space 
dependence. This scheme has recently been employed by 
Veenstra et al. 26]. TrBC allow one to limit the size of 
the space within which the numerical solution proceeds, 
without artificial (e.g. rigid-wall) boundary conditions. 
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FIG. 1: Demonstration of transparent boundary conditions. 
Solid line shows the integrated probability p{t) remaining in- 
side the computational domain; dashed line shows the sum 
of occupancy functions pv{t) over the set of basis functions 
{w„}l^i from sample A. 



In addition, modern computational power enabled us to 
perform more robust simulations which revealed aspects 
that have not been addressed before. 

For this work, an extension of the discrete TrBC as 
described in [25| was required to accommodate unequal 
saturation potentials on either end of the system (also 
applicable in case of a time-dependent potential in the 
inner region). Details of the finite element implementa- 
tion are described in appendix]^ To demonstrate perfect 
transmission through the domain border, we considered 
a Gaussian wavepacket with initial width of 50 nm, that 
was set free to slide down a linear potential ramp with 
slope F = 3 (Fig.[T]). As the wavepacket crosses the 
border, the total probability remaining in the domain 
falls smoothly to zero, showing that there is no reflection 
from the perfectly transparent boundary. 



B. Superlattice potential 

In order to avoid abrupt steps in the potential profile, 
we replaced the commonly used square barriers by an 
analytic form 
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X + a/2 X — a/2 
tanh — tanh 



(2) 



This is also a more realistic representation of an actual 
heterostructure potential [Tsj . 

We modeled a typical GaAs/Gai_:rAla;As heterostruc- 
ture (27} in the envelope function approximation. It has 
monolayer (ML) thickness of 0.283 nm and barrier height 
Vb = 790 X meV, x being the fraction of Al. The aver- 
age electron effective mass was set at m* = 0.071, to 
take account of non-parabolicity. The system has little 
sensitivity to the parameter a over the range 0.2 0.5 
nm. We chose a — 0.4 nm for all our samples, so that 
80% of the potential barrier height rises over two mono- 
layers. The characteristics of the potentials used in our 
simulations are laid out in Fig. [2] and Table HI 

Using TrBC assumes that the potential outside of the 
considered region is constant. Our investigations fo- 





FIG. 2: Band structure of the superlattice samples studied 
(from left to right: samples A, B, and C). 



Sample 


Vo, meV d, nm(ML) 


a, nm(ML) 


(T, nm 


A 


212 13.0 (46) 


3.1 (11) 


0.4 


B 


250 17.3 (61) 


2.5 (9) 


0.4 


C 


350 19.0 (67) 


2.5 (9) 


0.4 



TABLE I: Geometric parameters of the superlattice samples 
studied. Barrier height of 100 meV corresponds to 2;=0.13, of 
212 meV to a;=0.18, of 250 meV to a;=0.3 and of 350 meV to 
x=0.44 in the GaAs/Gai_a: ALAs structure. 



cussed on the dynamics of the wave packet inside the 
biased superlattice rather than those of the emitted part. 
The length of the computational domain, at least 40 cells, 
was sufficient for the simulations, since it allowed space 
for the packet to move while avoiding any contribution 
from surface states. Calculations over a wide range of 
bias with 40-cell and 120-cell domains showed no appre- 
ciable difference (less than 0.1%) in the output. Although 
the superlattices considered are ideal, one could intro- 
duce imperfections such as doping and barrier thickness 
fiuctuations, in order to study their dephasing effect on 
electronic coherence. 



C. Data analysis 

In this work, we will denote a Wannier-Stark (WS) 
level centered on the well with index k and belonging 
to miniband v as E^, the corresponding WS wavefunc- 
tion being Wl^{x) {u = 1,2,...); the well about which 
the initial wavepacket is centered is assigned index 0. 
By "resonance" (denoted as 71^^{X)) we will mean an 
anticrossing of energy levels and -E"^^" belonging to 
WSLz^ (the v^^' Wannier-Stark ladder, WSL) and WSL^, 
respectively, in sample X, where X can be A, B or C (so 
that E^ + nFd = where ij, fj. = 1,2,... and 

index n — 1,2, . . .). In terms of bias values, the term 
"resonance" will refer to a range of bias values that are 
close to the resonant bias {F = F-jzr,^^ or F„ for 7^^^,) and 
for which Rabi oscillations and/or RZT can be resolved. 
In case some of those indices are of little importance or 
have already been specified, they are omitted for brevity. 
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Throughout this work, the time unit is conveniently cho- 
sen to be the Bloch period at given bias: Tg = ^ and 
the length unit is the cell width d, unless stated other- 
wise. 

To visualize the interminiband dynamics of a 
wavepacket 'i'{x,t), we define the absolute occupancy 
functions 

pAt) = E (3) 

k 

which is the wavepacket intensity projected onto the i'*'* 
tight-binding (TB) miniband at time t (the index k runs 
over the cells within the computational domain) and the 
relative occupancy functions ^{t)= ^^^'^^^^j^ e [0,1]; 
j/=l,2,... 

Owing to the method of their construction, both the 
TB Wannier-Stark and Wannier states include the same 
harmonics, i.e. Bloch functions [28]. The main features 
of projection on minibands, such as resonant bias values 
and the Rabi oscillation period, using either set, were 
close to indistinguishable (with difference not exceeding 
1%) for the range of fields considered. Thus we adopted 
the simplification of using Wannier functions w^{x) as 
the projection basis. 

Generally speaking, the projection-on-minibands 
method does not apply at high fields where Wannier- 
Stark and miniband transport models [1] do not hold 
any more, and a sequential tunneling model has to be 
considered. Also, at the points of WSL energy level an- 
ticrossings, the considered TB WS functions [28|l can- 
not reliably describe resonant states (see e.g. [l^ llB])- 
However, at any bias the method gives a picture of the 
wavepacket 's distribution in energy (and hence between 
wells in real space), since w^lx) contain only harmonics 
with wavelengths A G [^d, ^^^^d]. For simplicity, we will 
use Wannier functions as the convenient orthogonal basis 
for miniband projection. 



III. RABI OSCILLATIONS: 
OVERVIEW 

A typical case of interminiband carrier dynamics is il- 
lustrated in Fig. [3] The time evolution of p^{t,F) is 
shown on a map plot comprising results from a large num- 
ber of time-dependent simulations over the range of bias 

= 0.1 . . .0.5 ^ with the step 1 a greater value 
of p{x, t) is shown in a lighter color. The large uniformly 
shaded areas correspond to exponential Zener decay of 
the wavepacket out of the superlattice potential. The 
darker vertical stripes originate from resonant Zener tun- 
neling where the wavepacket decays extremely quickly at 
energy level anticrossings. In the TB approximation, the 
n'^ anticrossing between WSL^ and WSL:/ occurswhen 
the condition Ei, — = F^nd = F^d is satisfied [29] . In 
Fig. [21 the period in -p of the series of stripes marked on 
top by shorter arrows implies interminiband separation 
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FIG, 3: (Color online) Absolute occupancy function map for 
the second miniband in sample A; '^(x,Q) — wiix). 



Fid — (223 ±5) meV and clearly corresponds to the anti- 
crossings associated with TZis, since E3 — E1 — 227,2 meV 
from our TB calculations, 

TZi3 do not show strong Rabi oscillations, since the 
energy gap following the third miniband, AE3, is only 
50 meV while an electron easily overcomes the intermini- 
band separation E2 — E1 = 89 meV, Therefore AE3 is too 
small to strongly bind an electron in the superlattice. On 
the other hand, AE2 = 73,02 meV suppresses tunneling 
to the third miniband at lower biases, and only 'R-\^^ are 
seen. 

There is also a set of prominent periodic spikes marked 
by longer arrows, corresponding to the group of reso- 
nances 'li\;2^- In contrast to T^is, the resonances TLi2 
do exhibit oscillations in p2{t) with time, corresponding 
to interminiband Rabi oscillations; they are wider and 
demonstrate a higher RZT rate for lower indices. The 
values =(6,9±0,2) ^ and Ft,^ =(2,33±0.05) S^Y 
are reasonably close to the anticrossing calculated in [lyj 
(7.2 and 2.4 respectively), even though the po- 
tential considered here was not exactly the square barrier 
one used in that work. 

Minor periodic changes in P2{t), creating a light hori- 
zontal mesh on the background with period Tb, are the 
signature of Bloch oscillations. For extremely high fields 
{F > 10 ^^^) the period and magnitude of these oscil- 
lations explodes (see e.g. the left edge of Fig. [TT|) . since 
AEi < F{d — a) and a transition to the next higher 
miniband can be made without timneling; then a theory 
different from interwcU hopping must be applied. 

Generally, it was found that energy anticrossings do 
not necessarily result in strong RZT for a strongly bound 
interacting WSL. Rabi oscillations appear to be the re- 
verse side of RZT in the WSL interaction: for strong 
RZT, overdampcd Rabi oscillations are seen (i.e. TZi^ in 
the above example), while persistent and strong Rabi 
oscillations correspond to weak RZT (e.g. 7^12). A 
quantitative relation between the two depending on the 
strength of the potential and resonance index remains 
an open question; up to now Rabi oscillations have been 
studied separately from RZT. 
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A. Resonance shapes 

The detailed structure of a resonance is shown in Fig. H] 
which is an enlargement of the resonance TZi2{A) from 
Fig. [3l Near a resonance, one sees Rabi oscillations as 
persistent oscillations of significant magnitude in p2(t) 
with period T-jz ~ 10 . . . 100 Tb (the corresponding fre- 
quency range is between the microwave and infrared re- 
gions). Both vertical and horizontal cross-sections of the 
plot demonstrate periodic oscillations of p2- 

A vertical cross-section of TZf2{A) at a near-resonant 
bias, shown at the bottom of Fig. [H explicitly demon- 
strates Rabi oscillations. Oscillations in pi{t) and P2{t) 
are shifted in phase by tt, which means that their cou- 
pling shows up as Rabi oscillations rather than RZT. In 
the example 72.^2 (^), the third miniband contributes lit- 
tle to total wavepacket norm. Its population could be 
due to two factors: (i) the non-TB WS functions of the 
second miniband involve harmonics from the third TB 
miniband, and (ii) the escaping part of the wavepacket 
passing through the third miniband on its way to the 
continuum. Since psit) oscillates in phase with P2{t), we 
conclude that the second factor is dominant in this ex- 
ample. 

When the asymptotic values of pi{t) and P2{t) in the 
limit t oo are not equal, as in the bottom panel of 
Fig. m we are dealing with a so-called asymmetric de- 
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FIG. 4: (Color online) Detailed view of 71^2 in sample A (top 
left). Period of Rabi oscillations vs. bias for TZl2{A) is fitted 
by means of a Lorentzian as in Eq. |4] (top right). Relative 
occupancy functions and wavepacket norm are shown for the 
bias F such that \^ — — I — 1.5 Ft^s ; pi is shown in solid 

"■12 

(red), p2 in dashed (green), pa in dotted (blue), p in chain- 
dotted (magenta) lines, and asymptotic occupancy values in 
horizontal lines (bottom). 
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FIG. 5: (Color online) Dynamics of the wavepacket having 
initial form 'l'(a;,0) — 'Wi{x) + wi{x + 3d) in reciprocal space 
at 7^?2(A). 



cay |30|]: the minibands 1 and 2 are coupled to the con- 
tinuum to a different extent. Right at a resonance bias 
value, the extremely strong interaction between the two 
minibands essentially merges them, and their asymptotic 
population values were both very close to 1/2. Away 
from the resonant bias, the coupling is weaker and the 
difference between the two values increases, eventually 
approaching unity. 

The dynamics of ^'(a::,i) in k-space demonstrates fea- 
tures similar to those of Bloch oscillations across a single 
miniband, i.e. most of the wavepacket steadily traverses 
the first and the second minibands as a whole with period 
2 Tb as shown in Fig. [5l The map plot consists of a se- 
ries of images of the total probability density p{x,t). We 
chose vE'(a;, 0) to be a linear combination of two Wannier 
functions in order to reduce the wavepacket uncertainty 
in k-space and to get a better resolution of fine details. 
From this perspective, an anticrossing can be thought of 
as a phenomenon of two adjacent WSL merging into a 
single broader one. 

For all resonance indices and samples examined, the 
Rabi oscillation period T-ji of an n*'' resonance clearly 
showed a Lorentzian-like dependence as anticipated from 
the nature of Rabi oscillations fsil : 



Tn{F) 



(4) 



which is shown in the top right panel of Fig. [H The 
data were obtained from a fit of P2{t) at TZ\2{A) over 
a long length of time considered [t — 50 . . . 100 T-ji). 
When speaking of half- width at half-maximum (HWHM, 
or Ftjii J, we will be referring to HWHM of the corre- 
sponding fit to T„(F) at -R^^. 

The maximum Rabi oscillation period observed at a 
resonant bias, T™"^, was found to grow exponentially 
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FIG. 6: Logarithmic fit of r,™"^ versus resonance index n for 
7^12(A) (left) and 7^?3(C) (riglit). 



with resonance index n in the following way: 



rpmax rpmax 

— ^1 



rpmax 
r^max 



(n-1) 



1,2,, 



(5) 



(see left part of Fig. [S]), which was ex- 
pected (32j since the transition matrix element 
\{W^{x)\W:p{x))\'^ (X e-l*^-™!. As the bias is re- 
duced, the horizontal tunneling channel to the next 
WSL lengthens, and it takes longer for the probability 
density to build up in the other miniband. We remark 
that perturbation theory predicts the dependence in 
Eq. (O to be linear [2l[, but that does not apply at the 
high bias considered. 



B. Pulsed output from the system 

Rabi oscillations of the carrier produce a periodic co- 
herent pulsed output with period T-ji similar to that of 
Bloch oscillations [3]. The data shown in Fig. [7] corre- 
spond to a record of jvl/endP = I'^ix end, t)\'^ values over 
time, at the endpoint x = Xend of the superlattice having 
a lower saturation potential. 

In the case of a relatively short Rabi oscillation period, 
the output consists of damped sinusoidal oscillations de- 
caying at the same rate as the wavepacket norm. For a 
long Rabi oscillation period accompanied by weak RZT, 
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FIG. 7: Pulsed output of a Rabi oscillating system at 
^12 (^) (top) and VJiziC) (bottom); chain-dotted line is the 
wavepacket norm. 
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FIG. 8: Dynamics of the electron dipole moment 
{'^{x,t)\x\'^{x,t)) calculated over the 0*'' (solid line) and 
the second (dashed line) cells at 7?.i2(>l). the wavepacket 
^(a:,0) = w-i(x) was taken at a near-resonant bias so that 
^ = 13.5 Tb (i = 0.296 -1^, i - ^ = 0.8 V^.J- posi- 

tion n is taken with respect to the center of the corresponding 
cell, in units of the cell width. 



as in the case of TZ\'^{C), the pulsed output is significantly 
distorted. Over the course of a single oscillation, the net 
effect of relatively small factors becomes significant; in 
the bottom section of Fig. [71 Rabi oscillations cause the 
drops in \^ end? with period Tn = 2050 Tb = 186 ps. 
Away from a resonance, the output pulses flatten out 
and on, a time scale longer than Tb, only a smooth ex- 
ponential decay in accordance with Zener theory would 
be observed. 

The form of the pulsed output is determined by the 
strength of Rabi oscillations, by their rate of magnitude 
decay, and by the strength of RZT. In principle, one could 
experimentally observe Rabi oscillations as well as estab- 
lish the relation between Rabi oscillations and RZT, by 
measuring the system's pulsed output. 



IV. WAVEPACKET SELF-INTERFERENCE 

A detailed consideration of the interference mecha- 
nism allows us to analyze resonance phenomena and draw 
some important conclusions without the need to calculate 
the system's spectrum. 

A strong correlation between subsequent well-to-well 
tunneling events arises due to their coherence. This 
correlation manifests itself in the interference between 
Bloch oscillations and intrawell oscillations (first resolved 
in [Toj ) with frequencies ujb and ujfn, that is constructive 
for near-resonant bias and destructive otherwise, with 
frequency detuning per single intrawell oscillation cycle 



Aw 



nuB 



(Efj. — El, 



Fnd 



(6) 



around a resonance with index n. If constructive, the 
interference leads to a gradual transition of the center of 
mass of probability density between the 0*'' and the n*^ 
cells and between the /z*'' and the i/*'' Brillouin zones in 
k-space. 
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At a resonant bias, the saturation value of occupancy 
functions of the two coupled minibands approached unity. 
Off resonant bias, the transition of probability density 
from one miniband to another was observed to reverse di- 
rection when a dephasing of vr had accumulated between 
the two resultant oscillations of wavepacket components 
residing in the two coupled minibands (Fig. [5]). One can 
clearly see intrawell oscillations (especially around t = 0) 
and Bloch oscillations, as well as the process of dephasing 
between Bloch and intrawell oscillations. If we call the 
wavepacket piece residing in the n*'' cell then with 
time, the phase shift between the peaks of the resultant 
oscillatory motion of ^I^o and ^'2 builds up, and reaches 
a net increment of tt at t = = 13.5 Tb- Remarkably, 
this interval corresponds to a maximum of 1^*2 (<)P and 
hence to a Rabi oscillation peak. 

Given these results, we conclude that interminiband 
Rabi oscillations are governed by the process of self- 
interference of a wavepacket subject to two intrinsic fre- 
quencies: ujb and w^^, and the amplitude of Rabi oscilla- 
tion at a given bias is determined by bias detuning from 
its resonant value. This again illustrates the coherence 
of multiwell tunneling, and in principle one can judge 
the coherence length of a superlattice by the number of 
observable resonances. 



A. Built-up state 

Let us consider the shape of the probability density 
that is building up, away from the wavepacket's origi- 
nal location, in the process of Rabi oscillation (built-up 
state). Without loss of generality, we will consider dy- 
namics of 4'(x, 0) = wi{x) at Tli2{A) and will denote by 
'i'Tz{x) a built-up state in the corresponding cell n. The 
shape of the probability density |\E'7?,(a;)|'^, turns out to be 
nearly the same for different values of the off-resonance 
field, and it also closely resembles the corresponding TB 
WS wavefunction, at least when |^'7j(a;)p is at its peak. 
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FIG. 9: Comparison of norms of a built-up state at F^2^(^) at 
resonant field for ^(a;, 0) = Wi{x) (dashed line) and '^'{x, 0) — 
1^1(0:) (dotted line) witli the corresponding Wannier-Stark 
state at TZi2 (A) (solid line line) in the left panel. In the right 
panel, the built-up state at a near-resonant bias is in dotted 
line. For reference, the biased superlattice potential is shown 
in double-dotted line. 



The latter comparison is made in the left panel of 
Fig. [51 Although in [l^ the corresponding computed 
resonant wavefunction had significant presence in the 



wells with indices 



7, our |4'7j(a;)| is not 



large there. A possible explanation is that the built-up 
state in the second miniband may include components 
from the third miniband as well. ^'-r,(x) appears to 
be only slightly more widespread than the TB W2{x), 
and on a linear plotting scale they are quite close as 
expected from [T^ for the moderate field considered. 
Note also that |^'K(a;)p shows slight asymmetry in k- 



at t 



— ^ — 13.5 Tb, there is more probability density 



space, in agreement with 15[: for example, in Fig. O 



present at negative values of k than at positive ones. This 
tells us that at moderate fields there is a great similarity 
between true and TB WS states. For comparison, the re- 
sults for ^'(a;,0) = Wi{x) are presented as well. Despite 
dissimilar shapes of the initial forms Wi{x) and wi{x), 
the built-up state looks almost the same for both. 

The built-up states at resonant and off-resonant val- 
ues of bias at 72.^2 are compared in the right panel of 
Fig. m The first one was taken at F such that the satu- 
ration population of the second miniband was p2 — 0.77 
{j; — jr — 0.5 r2). The second built-up state was taken 
at = F2 at the moment when the population of the sec- 
ond miniband reached the same value p2 = 0.77. Their 
close resemblance reveals that over the range of near- 
resonant fields, the build-up mechanism of |\E'K(a;)p is 
the same and produces a (rescaled) WS state of the sec- 
ond miniband at a resonant bias. 



B. Resonance condition 

In the process of Bloch oscillatory motion, the 
wavepacket tunnels out whenever it approaches the end 
of the Bloch oscillatory domain, producing a leaking out 
pulse. As a pulse propagates in space upon its escape, it 
scatters on the potential barriers, and some fraction of a 
pulse stays trapped in cells outside of the initial one. If 
the oscillations of the trapped part of a pulse happen to 
be in phase with those of a subsequent incoming pulse, 
the conditions for constructive interference are met and 
the probability density |5'7j(a::)p builds up in this well. 

For realistic fields, lub < i^^m (or Tb > and in 

order for oscillations of two subsequent pulses to be in 
phase in the n*'' well, the equality 



T^i/ 



UJB 



(7) 



must be satisfied: this is the condition for a resonance 
with index n. A gradual decrease in bias makes lob 
smaller, thus increasing n and shifting the location of 
the built-up state further down the potential ramp. Pro- 
vided that huifj^i, = En — and kujB = Fn d — ^d, 
for a resonance to occur between the i^*'' and the /j,*'' 



minibands, n 



In the TB approximation (E"^ 
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FIG. 10: Values ol Fiz^^nd plotted versus inverse bias for var- 
ious resonances with lowest indices (n = 1, 2, 3, . . . ascending 
from the left to the right) and compared with interminiband 
separation — Ev calculated from the tight-binding model. 
Filled circles mark the data corresponding to TZ\2{A), filled 
triangles to TZi2(C), empty circles to TZi^iC), empty rhombi 
to 72.23(C), empty triangles to TZi2{B). 



and are independent of F), that leads to a well-known 
result oc i [2^ . This relation is satisfied surprisingly 
well as indicated by results in the following subsection. 

From the above argument, it follows that the role of 
the interband jumping mechanism, as proposed in [30| 
for Rabi oscillation dephasing, must be minimal. Indeed, 
once tunneling of the carrier into a certain well (and re- 
flection back which doubles the phase increment of the 
ipnix, t)) puts it out of phase with the rest of the system, 
such well will not be largely populated due to destruc- 
tive interference. Mathematically, the discrete bias val- 
ues allowing constructive self-interference correspond to 
existence of poles of the system scattering matrix, which 
can be employed to successfully build a WS state for a 
multiband system [33 |. 



Resonant bias values 



that from the TB calculations {Efj_ — By), is less than 
2 meV or < 5%, for the resonances depicted. This is 
true even for resonances between non-ground minibands 
^23~^ (C) (this type of resonance has experimentally been 
observed in [39|). Thus even at high fields the equality 
F-ji-n^nd — Ey — Efj^ holds reasonably well, and the re- 
solved resonances occur periodically in y. However, in 
many cases this difference is seen to be much larger than 
the corresponding HWHM which cannot be explained by 
broadening of the transition line alone. Because F-jin jid 
is always larger than E"^ — Ey and this difference in- 
creases with bias, it must be a result of stronger coupling 
to the continuum at higher bias. This effect is similar to 
the energy level structure of a potential well becoming 
sparce as the well becomes shallower. 

Unlike other resonances, for the series 7^13 (C) the 
value of ndFn departs significantly from E^ — E2 for 
-p < 0.4 It has been verified that this deviation 

is not due to extreme narrowness of the first miniband in 
sample C. As follows from the series 72.23(C), the mutual 
alignment of levels belonging to WSL2 and WSL3 is al- 
most unchanged up to -p > 0.25; so it is the change in 
mutual arrangement of WSLl and WSL2 that is driving 
this deviation. 

The threshold value y « 0.4 corresponds to the 
potential drop per cell Fd « 47.4 meV which is close to 
the interminiband separation E2 — Ei = 44.6 meV. Ac- 
cording to the TB model, at |t < 0.4 the energy 
levels El and E2 would be found in the same cell. This 
is clearly impossible due to the Pauli exclusion principle, 
hence, the structure of separate Wannier-Stark ladders 
is completely destroyed at this point. At such bias, the 
Wi,{x) correspond only to the set of Fourier components 
with certain wavelengths rather than wavefunctions be- 
longing to certain minibands. 

Disruption of the WSL structure seems not to affect 
'^23(C), however: the initial wavepacket ^'(a;, 0) = ^2(2;) 
having only components with wavelength A G [f was 
noticed not to gain any significant components with 
A e [d, 2d] over time in our simulations, because of con- 
servation of energy. In general, the structure and prop- 
erties of 72.23 were found to be similar to those of TZu, at 
any bias. 



In order to check the validity of TB model calcula- 
tions for predicting resonant bias values at high fields, 
in Fig. (TU] the values of ndFn obtained from our simu- 
lations for the first few resonances in different samples 
are compared with the corresponding interminiband sep- 
arations calculated in the TB approximation, using a 
Kronig- Penney model. The data are plotted versus j:] 
x-error bars for all points refer to HWHM of a given res- 
onance r„ , whereas y-error bars show broadening of the 
energy levels and equal r„ rescaled in the same manner 
as i.e. ndTn- 

Excluding 72.13(C), the difference between the de facto 
relative position of the energy levels found as 7^7?,';,^ nd and 



V. RESONANCE ACROSS THREE MINIBANDS 

In a strong potential (sample C) having well-isolated 
minibands, we were able to resolve TZi'^'^{C); the loca- 
tion of several first resonances is indicated with longer 
arrows in Fig. [TT] To our knowledge, the phenomenon of 
Rabi oscillations across three minibands has been neither 
observed nor simulated before. 

Dynamics of Rabi oscillations at 72^3 shows that the 
wavepacket resides mostly in the first and the third mini- 
bands. In real space (upper panel of Fig.[T21 p^it) is seen 
as a part of p{x) whose shape has v humps per cell), there 
is little probability residing in wells -1 and -2 at all times. 
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In reciprocal space (lower panel), the second Brillouin 
zone consistently remains underoccupied relative to the 
first and third zones; some traces of p{k,t) in the fourth 
and the fifth Brillouin zones correspond to RZT, and the 
finer background oscillations are caused by Bloch oscilla- 
tions in individual minibands. Thus, the carrier mostly 
bypasses the second miniband since the resonance con- 
dition for TZi2 is not well satisfied, and tunnels directly 
into the third miniband. 

It was found that TZi2 and 7^13 have many features 
in common. Namely, the resonance shape was similar, 
the oscillations of pi{t) were close to sinusoidal, Eq. (|4]) 
was satisfied, and in k-space most part of a wavepacket 
^{x,0) = wi{x) traversed the first three Brillouin zones 
of the TB model (fc G [-^, ^]) as a whole. However, 
individual intrawell oscillations were not well pronounced 
at TZi3 because of the strong potential barriers in sam- 
ple C and the fact that intrawell oscillations between the 
pairs of coupled minibands l<->2 and l<->3 are of compa- 
rable magnitude and strongly interfere with each other, 
as was noticed from the irregular shape of the result- 
ing oscillations of a wavepacket within the 0*'' cell. Also, 
within the 0*'' cell, the center of the resultant oscillations 
of the probability density was shifted down the potential 
ramp from the center of the cell, under the influence of 
strong bias which creates asymmetry about the cell cen- 
ter. 
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A. Role of sandwiched miniband 

For a resonance across three minibands, an exponential 
fit based on the two-miniband approximation is not a 
good fit to T™°^ (n) any more (right panel of Fig. [5]). The 
reason lies in the involvement of the second miniband 
which is "sandwiched" between minibands 1 and 3, in 
the interminiband tunneling process. Despite its small 
average population value, the second miniband must be 
taken into account at any bias as will be shown below. 

The role of the second miniband in carrier transfer be- 
tween minibands 1 and 3 must depend on resonance in- 



f-f^ ttnS TP J 'JJ^ 




II. I 0,2 t),.1 (1.4 11,5 (1.6 (17 (l,K (t,') I 
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FIG. 11: (Color online) Absolute occupancy function map for 
the first miniband in sample C; 4'(a;,0) — Wi{x). 



FIG. 12: (Color online) Dynamics of the wavepacket 
^{x,0) — 'Wi{x) in real (top) and reciprocal (bottom) space 
at TZlg in sample C. 



dex, n. Otherwise, {p2{t)) should be proportional to the 
interminiband tunneling rate of '^{x,t), or inversely pro- 
portional to the period of Rabi oscillations, and the latter 
would be described by Eq. ([5]). Hence, this would imply 
{P2) oc e~". 

Values of (^2) obtained are shown in Fig. [131 which 
displays dynamics of the relative occupancy function at 
7^13 (C). It is remarkable that, despite the strong cou- 
pling to the third miniband, the RZT rate is small for 
the resonance indices n > 2 thus allowing for long-lasting 
Rabi oscillations. For lower resonance indices, the popu- 
lation of the second miniband reaches significant level at 
certain times. As the resonance index ascends, the aver- 
age population of the second miniband at 7?,"3 decreases 
(Fig. [131); its population is largely caused by Rabi oscilla- 
tions between minibands 1 and 3. As p2{t) 0, most of 
the wavepacket undergoes sinusoidal oscillations between 
minibands 1 and 3, i.e. for lower biases the interminiband 
dynamics resembles a two-miniband model (which would 
assume P2it) = 0). The values of (^2) are shown sepa- 
rately in Fig. [TH and appear to significantly deviate from 
exponential dependence at TZf^. 

In the process of Rabi oscillations at TZis, the popula- 
tion level of the sandwiched second miniband is higher, 
the better is the match in energy between the states from 
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FIG. 13: (Color online) Occupancy functions and wavepacket 
norm at TZi-j{C) for n = 1 (top) through 5 (bottom), pi in 
solid (red), p2 in dashed (green), ps in dotted (blue)and p 
in chain-dotted (magenta) lines; the broken black line shows 

(P2>. 



WSL2 and the initial carrier energy E^. The reason lies 
in the tunnchng mechanism involved: in the process of 
the multiwell tunnehng between minibands 1 and 3, an 
electron's energy after a certain number of interwell hops 
can be close to i?2- This proximity greatly enhances pop- 
ulation of the sandwiched miniband by providing avail- 
able density of states. Conversely, a large population 
level of a sandwiched miniband favors transition between 
the first and the third minibands. When the transi- 
tion element a^il becomes comparable to (here 



{w]l; {x)\x\w^ (x))) , indirect tunneling through the 
sandwiched miniband becomes significant. Therefore the 
degree of participation of the second miniband and its ef- 
fect on 7^13 is determined by alignment of energy levels in 
WSLl, WSL2 and WSL3 and is dissimilar for resonances 
with different indices. The details of the influence of this 
alignment on the electron dynamics will be given in the 
next subsection IV Bl 

The correlation between coupling strength to the sand- 
wiched miniband and the rate of interminiband tran- 
sition, or Rabi oscillation frequency, is clearly demon- 
strated in Fig. [131 for TZf^, both ip- and (^2) are smaller 
than expected from the exponential fit. For index n — 4, 
the tunneling channel between the levels and £'| con- 
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FIG. 14: Rabi oscillation frequency (left scale, filled circles) 
and (P2) (right scale, empty circles) at 7?.i3(C) fitted with 
exponential curves (chain-dotted lines) on a logarithmic scale. 



sists of four interwell hops, and the intermediate energy 
levels available to the carrier are E''. In this case, the 



^2 

alignment of i?2'^ turns out to be particularly unfavor- 



able for them to act as a strong transition channel. At 
F F4 = 1.245 E^ and E^ are equally remote 



1.245 

mcV ' 

from the carrier's initial energy: {E^ — I ~ l-^i ' 



« 20 meV (whereas in the other cases, when n ^ A 
and n > 1, one of the levels E2 is closer to than 7 
meV). Such an alignment can also be anticipated from 
Fig. [HI where TZ\^ lies in the middle between TZ\2 and 
TZ\2 and hence is particularly isolated from any of the 

-ryk 
'^12- 

A sandwiched miniband affects not only Rabi oscilla- 
tions, but also RZT. In Fig. [31 we can see that despite a 
weaker bias, RZT at Ti,\^ is stronger than for the preced- 
ing T^fg for the same reason, namely proximity to Ti-\2- 
There, it also makes Rabi oscillations corresponding to 
TZi2 vanish at the given field since the barrier through 
which the carrier tunnels to the third miniband and then 
to the continuum is significantly reduced by coupling be- 
tween WSL2 and WSL3. In other words, closely situated 
resonances mutually influence each other and can be re- 
ferred to as coupled resonances. 

Despite the sometimes small value of its average pop- 
ulation, the second miniband has a significant effect on 
resultant interminiband dynamics. Thus to calculate pa- 
rameters of a resonance across three minibands, it is nec- 
essary to take at least these three minibands into consid- 
eration. 



B. Coupled resonances 

In order to understand the interaction between differ- 
ent resonances better, let us consider an example where 
strength of interminiband coupling for different reso- 
nances is comparable and the resonances strongly inter- 
fere with each other. The region corresponding to the 
coupling of two resonances Ti-2z and TZ\2 in sample B is 
shown in Fig. 1151 a map plot of p2 for 4* (a;, 0) = 'W2{x). 
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71^2 is situated on the left and 7?.23 on the right, the 
former being wider due to its lower resonance index. 

Detail of the carrier dynamics for field ^=0.646 
where contributions from T^fj ^^'^ ^23 nearly equal, 
are shown in the bottom section of Fig. [151 Since all three 
minibands are interacting at the same time, the resultant 
dipole dynamics looks somewhat similar to a superposi- 
tion of the three corresponding Rabi oscillations for the 
individual resonances {TZi2, 7^23 ^13); ^^'^ time 
evolution of pi,2,3(i) exhibits beats. Therefore for cou- 
pled resonances, the resultant Rabi oscillations are an 
interference product of not only Bloch and intrawell os- 
cillations, but also of Rabi oscillations corresponding to 
the other pairs of strongly coupled minibands in the sam- 
ple. This interference allows the carrier to use the "extra" 
density of states to enhance tunneling. A map plot of the 
wavefunction in the coupled resonances zone (Fig. I16p 
demonstrates that large miniband populations facilitate 
the build-up of each another. For example, p2(t) is rel- 
atively low at < = 35 Tb and t = 70 Tb, i.e. at the 
moments when P3{t) is almost zero. Likewise when p2(t) 
is relatively large, psit) becomes larger as well (e.g. at 
t = 20 Tb, 50 Tb and 90 Tb). 
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FIG. 15: (Color online) Detailed view of coupled resonances 
(7?.i2,7?.23 and TZl^) in sample B for ^'(j:, 0) = W2{x). Relative 
occupancy functions dynamics at bias -i =0.646 is shown 
in the bottom panel; pi is shown in solid (red), p2 in dashed 
(green), p3 in dotted (blue), and p in chain-dotted (magenta) 
lines. 
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FIG. 16: (Color online) Dynamics of the wavepacket 
$(a;,0) — wi{x) in direct space in the coupled resonances 
zone in sample B. 

Obviously the two-miniband model cannot reproduce 
such behavior. As was recently found [34j, Rabi oscilla- 
tions in a three-level system are not merely a summation 
of Rabi oscillations between the three separate pairs of 
minibands, they rather are a coherent superposition. In 
fact, the temporary reduction of Rabi oscillation ampli- 
tude observed in work was due to Rabi oscillation 
revival, i.e. the beats caused by superposition of Rabi os- 
cillations between several lowest minibands, rather than 
to Rabi oscillation dephasing. 

VI. CONCLUSION 

In conclusion, our time-dependent simulations of elec- 
tron dynamics in a biased superlattice provide an 
overview of carrier behavior over a large range of bias. 
We were able to identify energy level anticrossings as 
resonances displaying Rabi oscillations and/or RZT and 
to study the structure of these resonances in detail; our 
findings could be experimentally studied by observing the 
system's pulsed output. 

It has been shown that Rabi oscillations result from 
constructive interference between Bloch and intrawell os- 
cillations. A distinct interminiband resonance occurs 
whenever the initial wavepacket is capable of building 
up a Wannier-Stark state of significant magnitude away 
from its initial location, through the coherent process of 
self- interference. 

We have also reported and analyzed Rabi oscillations 
across three minibands. The important role of a sand- 
wiched miniband in the transitions across more than one 
miniband has been studied: energy-wise it provides avail- 
able density of states in the intermediate tunneling re- 
gion, and dynamically it invokes additional Rabi oscilla- 
tions interfering constructively with the principal inter- 
miniband transition. 

At a resonance across three minibands and in the cou- 
pled resonances region, the resultant carrier dynamics 
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has significant contribution from all three minibands. 
Hence these cases cannot be described within the com- 
monly used two-miniband model, and require taking at 
least three minibands into consideration. 

Interminiband separations calculated in the tight- 
binding approximation can be used to predict resonant 
values of bias in the high-field regime with reasonable 
accuracy for a variety of potentials, including resonances 
between minibands 2 and 3 and those across three mini- 
bands (e.g. TZi3, with the only restriction Fd < E2—E1). 

Although the present work focuses on an ideal zero- 
temperature superlattice under uniform constant elec- 
tric field, the numerical methods utilized are easily ca- 
pable of handling time-dependent irregular potentials. 
With minor variations, the numerical techniques can 
be applied to study carrier dynamics in other systems 
(e.g. double quantum dots [3, [s^ and irregularly 
shaped potentials 36]) and even areas of physics, such 
as photonics ,37. , ,38, . Isol . |40| and cold atom optical 
traps [41I, m, |43[ |44| : these are reserved for the future 
investigations. 
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APPENDIX A: DISCRETE TRANSPARENT 
BOUNDARY CONDITIONS 

Our numerical method was based on TrBC imple- 
mented with Numerov and Crank-Nicholson methods for 
space and time respectively, with the cumulative preci- 
sion 0{{6x)^) in space and 0[{5tY) in time. We ex- 
tended the discrete TrBC as described in [2^ to the case 
of unequal saturation potentials on opposite sides. For 
the most part, we keep the same notation as in (25j . 

Let us express Eq. ^ in terms of finite differences. We 
will start from the finite difference form of the system's 
propagator which translates the time by 5t 



«'(a;,t + A) 



-iH(t)A 



with A = ^ and H{t) being the (possibly time- 
dependent) Hamiltonian of the system. Cayley's approx- 
imation preserves unitarity and is exact to second order: 



-,-iHA 



1 



iHA 



1 + ^iHA 



A little algebra leads to the expression 
2 mm 



V{x,t) 



2i\ 

~k} 



(*(x,f + A) + 
^ix,t) + 0(A3) 



with the notation V{x,t) = Vsl{x) + xF. It serves as 
a discretized in time version of Eq. ([T]) and as a start- 
ing point for discretization in space using the Numerov 
method. 

In order to build the solution 5'(a;,<) to Eq. ([T]), we 
will use a uniform space grid consisting of J points and 
defined as Xj = j Sx, (j = 0, 1, . . . , J — 1) and uniform 
time grid i'-'^-' = k St, {k — 0,1, . . .). The subscript j will 
refer to the inner points Xj , and e to the end points xq or 
as appropriate; the superscript (fc) will denote the 
time instant t'^'^K Thus '^{xj,t^'''>) = '^f ^ and V{xn) = 

To apply TrBC in the Numerov approximation for the 
ID case, the necessary conditions are: V{x < xq) = 
V{x — xo) = Vo, and V{x > xj^i) = V{x = xj^i) = 
Vj-i, as well as = at both ends. Having started 
at i = 0, we proceed as follows to construct ^(x, i^'^''"^)) 
from ^{x,t'^''^): 

(i) Calculate the time- independent coefficients over the 
space grid: 



9j 



2 mm* 2i 

~ir~ ^ A ^ 



1 



{5xf 



12 



■9j 



2 + 



< j < J - 1 



< j < J - 1 
< j < J - 1 



(ii) Calculate the time-independent border coefficients 
(e = or J — 1 separately): 



1 



2 de 



1 



2 i {dx)'^mm* 1 



3St 

'at - 1^ 



Ae = 
Pe = 



arg 

1 - la, 



\l-ai\ 
de{ae - ae) 
d*{al - ae) 



(iii) Construct the polynomials for the next time step: 
p(.) ^ 



1 

Ae 
2k+l 



fe+l "^e^ e 

1 p~i(k- 
2k+l ^ 



4 Pi'^-l) ^ArPe 



(fe-2) 



k+1 



k=-l 
k = 

k>0 



(fc) 



k^O 

t-'^) k>0 
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(iv) Calculate the time-dependent coefficients over the space grid: 
I 



ISTT + i^r 0<,<J-1 



[k+l 
9,7-2 



1 - aj_iej_2 



J = ^ - 1 
< j < J 



(v) Based on the above, construct wavefunction for the 
next instant of time: 



(fe+l) _ '^3 



and Pe'^'^ being Legendre polynomials of the /c"* order 
having A,, at the time instant t = t^'^^ as an argument. 



with the notation 



c 

A 



_ 8i (Sx)'^ mm* 

" St 

_ 2i ((5x)^ mm* 



3 6t 



St 



If the potential considered is time-dependent in the 
inner region, the coefficients from step (ii) will change 
with time, and we will simply have to recalculate them 
for every new instant of time. In the formula language, 
this means substitution everywhere above V{xj) 

V{xj,t^'''>) = V- '^\ which results in time-dependent coef- 



ficients gj — > 5^ , d 



etc. 
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